Finite-temperature phase diagram of nonmagnetic impurities in high-temperature 
superconductors using a d = 3 t J model with quenched disorder 
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We study a quenched disordered d — 3 tj Hamiltonian with static vacancies as a model of 
nonmagnetic impurities in high-T c materials. Using a renormalization-group approach, we calculate 
the evolution of the finite-temperature phase diagram with impurity concentration p, and find several 
features with close experimental parallels: away from half-filling we see the rapid destruction of a 
spin-singlet phase (analogous to the superconducting phase in cuprates) which is eliminated for 
p > 0.05; in the same region for these dilute impurity concentrations we observe an enhancement 
of antiferromagnetism. The antiferromagnetic phase near half-filling is robust against impurity 
addition, and disappears only for p > 0.40. 
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The electronic properties and phase diagram of high- 
Tc materials are particularly sensititive to impurities — 
substitution of 3d transition elements (Zn, Ni, Co, Fe), 
or other metals (Al, Ga), for the Cu atoms of the Cu02 
planes The interplay between disorder, strong anti- 
ferromagnetic correlations in the parent compound, and 
doped charge carriers, offers a window onto the nature 
of both the superconducting phase and the normal state 
above T c . Doping by nonmagnetic (S — 0) Zn ions pro- 
vides one representative example: the most pronounced 
effect is the rapid destruction of the superconducting 
phase [HQ; m YBCO the transition temperature is re- 
duced at a rate of ~ 15K/at.% of impurities, so that 
it takes Zn concentrations of only about 6% to entirely 
eliminate superconductivity |2| . This is in contrast to the 
antiferromagnetic phase at half-filling, which requires a 
far larger Zn concentration (about 40% in LSCO [3[) to 
completely suppress. The effects in the metallic region 
above T c are equally surprising: nuclear magnetic res- 
onance experiments have found that Zn atoms induce 
local magnetic moments at nearest- neighbor Cu sites [3] , 
and enhance antiferromagnetic correlations for several 
lattice spacings around the impurity 0, Q. In lightly 
hole-doped LSCO there have been observations of an ini- 
tial increase in the Neel temperature with Zn addition, 
and even impurity-induced reappearance of long-range 
antiferromagnetic order 0, [1] . 

In this work we model the effects of nonmagnetic im- 
purities in high-T c materials through a d = 3 t J Hamilto- 
nian with quenched disorder in the form of static vacan- 
cies. Through a renormalization-group (RG) approach 
we obtain the evolution of the global temperature ver- 
sus chemical potential phase diagram with disorder. Our 
results capture in a single microscopic model some of 
the major qualitative features of impurity-doping in real 
materials: the rapid suppression of a spin-singlet phase, 
analogous to the superconducting phase in cuprates; the 
gradual reduction of the antiferromagnetic phase near 
half-filling; and the enhancement of antiferromagnetism 
away from half-filling for small impurity concentrations. 



We consider the quenched disordered tJ model 
on a d-dimensional hypercubic lattice, —j3H = 

T,(ij){-P H o(i,j)} + Y,il*T Vn ii where -P H o(iJ) = 
~t E ff ( c L c j* + c] a c ia ) ~ J (Si • Sj - UiUj /4)+n(n l + rij) 
is the standard t J model pair Hamiltonian. The static 
impurities at each site i occur with probability p via 
/i™ p = — oo and do not occur with probability 1 — p 

imp p. 

via n i =0. 

To formulate an RG transformation for this system, 
we use the d = 1 Suzuki- Takano decimation [9l-[l9|. gen- 
eralized to d > 1 through the Migdal-Kadanoff method 
[20l [2l1 |. This technique, adapted for quenched random 
bond disorder, has recently elucidated the phase dia- 
grams of the quantum Heisenberg spin-glass in various 
spatial dimensions [18j . In our case the rescaling for the 
d = 1 system (with sites z = 1,2,3,...) is: 
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where the traces and sums are over even- or odd- 
numbered sites i, and —/3'H' is the renormalized Hamil- 
tonian. Anticommutation rules are correctly accounted 
for within segments of three consecutive sites, at all suc- 
cessive length scales as the RG transformation is iterated. 

The algebraic content of the RG transformation is con- 
tained in the second and third lines of Eq. ([1]), yielding 
the renormalized pair Hamiltonian — ft' H' a {i' ', f) through 
the relation: exp(- /3' H' (i' , /)) = Tr fc exp(-(3H (i' , k) + 
/ij™ rifc — /3Ho(k,j')). Under the transformation the 
original system is mapped onto one with a more 
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FIG. 1: Hierarchical lattice on which the d — 3, b — 2 Migdal- 
Kadanoff recursion relations are exact. 



general form of the pair Hamiltonian, — /3Ho(i,j) 



where the interaction con- 
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distributed with a joint quenched probability distribu- 
tion ■p(Kjj). This generalized form of the Hamilto- 
nian remains closed under further RG transformations. 
Through the relation above we can write the interac- 
tion constants K^, of the renormalized pair Hamil- 
tonian — f3' H' (i' , f) as a function of the interaction 
constants K^/^ and Kfc,j' of two consecutive nearest- 
neighbor pairs in the unrenormalized system, K^, •, = 
R(Ki'fe, Kfcj/). This function R comes in two 
varieties, depending on whether or not there is 
an impurity at site k, which we shall denote as 
Ro and Rj mp respectively. Starting with a sys- 
tem with quenched probability distribution 'P(Kjj), 
the distribution 7"(K^, -, ) of the renormalized sys- 
tem is given by the decimation convolution (22j: 
V'{K' Vjl ) = JdK i , k dK. kj ,V{K i , k )V(K kj ,)\p5(K^ JI - 
R^ p (K i , fe ,_K fe/ ))+(l-p)5(K^, r Ro(K ia ,K^))]. The 
initial condition for the RG flow is the distribution corre- 
sponding to the original system, 'Po(Kjj) = S(Kij — Kq), 
where K = {t, J, - J/4, fi, 0}. 

The RG transformation is extended to d > 1 through 
the Migdal-Kadanoff 0, Hl| procedure. While approxi- 
mate for hypercubic lattices, the recursion relations gen- 
erated by this procedure are exact on hierarchical lattices 
[23U2H , and we shall use this correspondence to describe 
the RG transformation for the case d — 3, with length 
rescaling factor 6 = 2. The associated hierarchical lattice 
is shown in Fig. 1. Its construction proceeds by taking 
each bond in the lattice, replacing it by the connected 
cluster of bonds in the middle of Fig. 1, and repeating 
this step an infinite number of times. The RG trans- 
formation consists of reversing this construction process, 
by taking every such cluster of bonds, decimating over 
the degrees of freedom at the four inner sites of the clus- 
ter, yielding a renormalized interaction between the two 
edge sites of the cluster. Denoting these edge sites as i', 
j' , and the four inner sites as k\, . . . , fej, this decimation 
can be expressed as Kj/j/ = Xm=i ■^(■K-i'fcn^fcni')' ^ ust 
as in the d — 1 case, this decimation will give, after a 
single RG transformation, a system with a nonuniform 
quenched distribution of interaction constants. We can 
calculate the quenched probability distribution V'(K.iij>) 
of the renormalized system through a series of pairwise 
convolutions, consisting of the decimation convolution 
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FIG. 2: Pure system (p = 0) phase diagram of the isotropic 

d = 3tJ model [UG3I for J/ t = °- 444: ( a ) in terms of chem- 
ical potential jj,/J vs. temperature 1/t; (b) electron density 
(rii) vs. temperature 1/t. Panels (c) and (d) show the analo- 
gous phase diagrams for the uniaxially anisotropic case [l7| ]. 
with t z /t xy = 0.3, J z /J xy = 0.09, J xv /t xy = 0.444. In both 
cases antiferromagnetic (AF), dense disordered (D), dilute 
disordered (d), and r phases are seen. The solid lines rep- 
resent second-order phase transitions, while the dotted lines 
are first-order phase transitions (with the unmarked areas in- 
side corresponding to coexistence regions of the two phases at 
either side). Dashed lines are not phase transitions, but disor- 
der lines between the dilute disordered and dense disordered 
phases. 



defined above for interactions in series, and a "bond- 
moving" convolution for interactions in parallel, using 
the function Rb m (K^4,Ks) = + Kg. In order to 
numerically implement the convolution, the probability 
distributions are represented by histograms, where each 
histogram is a set of interaction constants (t, J, V, fi, v) 
and an associated probability. Since the number of his- 
tograms that constitute the probability distribution in- 
creases rapidly with each RG iteration, a binning pro- 
cedure is used [26|. Furthermore since evaluation of the 
R functions is computationally expensive, and most of 
the weight of the probability distributions is carried by a 
fraction of the histograms, we have added an additional 
step before the decimation convolution to increase effi- 
ciency: the histograms with the 100 largest probabilities 
are left unchanged, while the others are collapsed into a 
single histogram in a way that preserves the average and 
standard deviation of the quenched distribution. Thus 
we evaluate 10 4 local decimations at each RG transfor- 
mation. 

All thermodynamic properties of the system, in partic- 
ular the finite-temperature phase diagram, can be deter- 
mined from analyzing the RG flows. In the pure (p = 0) 
case, the transformation described above reduces to the 
recursion relations derived for the d = 3 tJ model in ear- 
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FIG. 3: On-site and nearest-neighbor singlet and triplet pair- 
pair correlations for the d = 3 tj model, with p = 0, J/t = 
0.444. In (a) and (b) the correlations are plotted as a function 
of chemical potential fi/J at constant temperature 1/t = 0.10. 
In (c) they are plotted as a function of temperature 1/t, at 
the constant electron density (n/) = 0.67. The corresponding 
phases are indicated near the top of each plot, with solid 
and dotted vertical lines marking second-order and first-order 
phase boundaries respectively. 



lier studies [Tl|, [l2| , and yields the phase diagram shown 
in Fig. [H[a,b) for J/t = 0.444. Here we summarize the 
observed phases (for details, consult [HI, [Hi) : near na, lf~ 
filling (/i/J — > oo, (rii) — > 1), there is a transition with 
decreasing temperature from a densely-filled disordered 
phase (D) to long-range antiferromagnetic order (AF). 
This AF phase persists away from half-filling down to 
l±/J~ 1.6, or 5% hole doping. For very large hole dop- 
ings (> 37%) we go over into a dilute disordered phase 
(d), with narrow first-order coexistence regions between 
the d and D phases. At intermediate hole dopings of 
33-37% a novel phase (r) is found at low temperatures, 
flanked by an intricate lamellar structure of AF islands. 

The r phase is characterized by the formation of 
nearest-neighbor spin-singlet pairs, as can be under- 
stood from correlation functions calculated using the 
RG flows. Let us define a singlet pair-pair correla- 
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utive sites i, j, k in the lattice, Fig. [3] shows the on-site 
correlations T^J p T^g, and nearest -neighbor correla- 
tions T^ e k , T*pj k . In Fig.GUa) and (b), we see a constant 
temperature slice at 1/t = 0.10 as fi/J is varied. There 
is a broad region of chemical potentials away from half- 
filling, centered at the r phase, where both the on-site 
and nearest-neighbor singlet correlations are strong, in 
contrast to the triplet correlations, which are suppressed 
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FIG. 4: Calculated phase diagrams of the d = 3 tj model, 
with J/t — 0.444, for various values of the impurity concen- 
tration p, plotted in terms of temperature 1/t vs. chemical 
potential fi/J. The phases depicted in the figures are: dilute 
disordered (d), dense disordered (D), antiferromagnetic (AF), 
and r. The inset shows AF transition temperatures for the 
near-half-filled system (p./J = 100) as a function of p. 



in the same region. We see similar behavior in Fig. EHc), 
where the correlations are plotted as a function of tem- 
perature 1/t at a constant electron density (rij) = 0.67. 
As we decrease the temperature, approaching the transi- 
tion into the r phase, there is a significant increase in the 
singlet correlations and rapid decay of the triplet correla- 
tions. Spin-singlet liquids, i.e., the hole-doped resonating 
valence bond (RVB) state, have featured prominently in 
theories of high-T c superconductivity (for a review see 
Ref. |27[). As we shall see below, the behavior of the r 
phase under impurity-doping is analogous to that of the 
superconducting phase in high-T c materials. 

Though in this study we focus on the isotropic d = 3 
model, there is evidence that general features of the 
p = phase diagram discussed above persist even in 
the case of spatial anisotropy, with uniform interactions 
(Jxy, Jxy) along the xy planes and weaker interactions 
(t z , J z ) along the z direction. Through a similar RG 
approach, using the more complicated hierarchical lat- 
tice associated with a uniaxially anisotropic cubic lat- 
tice [28| , it was found in particular that the r phase con- 
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tinues to exist in roughly the same doping range even 
for weak interplanar coupling, though as expected the 
transition temperatures steadily decrease as the coupling 
is reduced [17|. A representative phase diagram, with 
tz/txy = 0.3, J z /J xy = 0.09, J xy /t xy = 0.444, is shown 
in Fig. [2][c,d). Thus the r phase may be relevant even 
in the strongly anisotropic regime important for high-T c 
materials, which are characterized by weakly interacting 
CuC>2 planes. 

In Fig. 2] we show the evolution of our calculated phase 
diagram with increasing impurity concentration p. The 
r phase is rapidly suppressed for p = 0.01 through 0.04 
[Fig. H[a)-(d)], and is no longer present by p — 0.05. 
The rate at which the r phase disappears is comparable 
to the reduction of T c with nonmagnetic impurities in 
cuprates, where typically concentrations w 2 — 6% (de- 
pending on dopant) are enough to eliminate supercon- 
ductivity (HQ. As the area of the r phase recedes for 
these small impurity concentrations, the region it for- 
merly occupied is replaced by a complex lamellar struc- 
ture of the AF phase. We can understand this enhance- 
ment of antiferromagnetism through an RVB-like picture 
of the t phase [2!| : in the pure case the nearest- neighbor 
singlets resonate in all possible arrangements along the 
bonds; when an impurity is added some of these arrange- 
ments are "pruned", because the bonds adjacent to the 
impurity can no longer accommodate singlets. This in- 
hibition of singlet fluctuations leads to enhanced anti- 
ferromagnetic correlations around the vacancy. Such lo- 
cal AF enhancement near dilute nonmagnetic impurities 
has been observed through NMR and NQR studies on 
Zn-doped YBCO @, and supported theoretically by 
finite-cluster studies of the d = 2 Heisenberg [3l[ and t J 
HHHl] models. More dramatically, in lightly hole-doped 
La2- x Sr s Cui_ 2; Zn z 04 (with x = 0.017) the Neel tem- 
perature actually increases with the addition of Zn up to 
z = 0.05, before turning downwards again at higher z Q- 
A similar, though smaller, effect has been found even at 
larger hole dopings of a; = 0.115 and 0.13, with the Tjy 
increasing up to z = 0.0075 8]. In the case of x = 0.13, 
there is even no long-range antiferromagnetic order for 
the Cu spins in the Zn-free compound; it appears for 
z > 0.0025. This reappearance of long-range AF or- 
der upon addition of impurities, at small hole-dopings 
away from half-filling where it does not exist in the pure 



case, was replicated in the d = 2 t J model using a self- 
consistent diagrammatic approach [34j , and in the d = 2 
Hubbard model with the dynamical cluster approxima- 
tion [35[. Thus the enhancement of the AF phase away 
from half-filling, which wc find at small impurity con- 
centrations, is consistent with previous experimental and 
theoretical indications. 

On the other hand for larger concentrations of impu- 
rities, the dilution of the spins in the lattice becomes the 
dominant effect, and eventually all long-range magnetic 
order is destroyed in the system. We see this in Fig. 0|e)- 
(h), showing phase diagrams for p = 0.10 through 0.40, 
and in the inset which plots the AF transition temper- 
ature as a function of p near half-filling (/i/J = 100). 
In contrast to the r phase, the AF phase around half- 
filling is robust against impurity addition, and only dis- 
appears for p > 0.40. Qualitatively similar behavior has 
been seen in the half-filled compound La2Cui_ z Zn z 04, 
where Zn concentrations of z s=s 0.4 are required to reduce 
the Neel temperature to zero 3j , much larger than those 
needed to eliminate superconductivity in the hole-doped 
material. 

To summarize, we have applied an RG approach to 
the quenched disordered d — 3 tJ model, and found the 
evolution of the phase diagram as a function of impu- 
rity concentration. The spin-singlet phase away from 
half-filling is quickly destroyed through the addition of 
small quantities of static vacancies, while antiferromag- 
netism in the same region is enhanced. The antiferro- 
magnetic phase near half-filling is less sensitive to impu- 
rity addition, and completely disappears only at larger 
concentrations. These results all have close parallels in 
experimental results from cuprates. The RG method de- 
scribed here for dealing with quenched disorder in the t J 
Hamiltonian could be generalized to more complex sys- 
tems: for example the disordered Hubbard model, where 
the double-occupation of sites is allowed through a finite 
electron-electron repulsion. The role of electron correla- 
tions and disorder in this system has led to interesting 
phase diagram predictions |36l - l38l ]. which could be fur- 
ther explored with RG techniques. 
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